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1. Introduction and Summary 

One of the most challenging problem is to find consistent formulation of massive gravity. 
The first attempt for construction of this theory is dated to the year 1939 when Fierz and 
Pauli formulated its version of linear massive gravity [0 . However it is very non-trivial 
task to find a consistent non-linear generalization of given theory and it remains as an 
intriguing theoretical problem. It is also important to stress that recent discovery of dark 
energy and associated cosmological constant problem has prompted investigations in the 
long distance modifications of general relativity, for review, see Q. 

It is natural to ask the question whether it is possible to construct theory of massive 
gravity where one of the constraint equation and associated secondary constraint eliminates 
the propagating scalar mode. It is remarkable that linear Fierz-Pauli theory does not suffer 
from the presence of such a ghost. On the other hand it was shown by Boulware and Deser 
H that ghosts generically reappear at the non-linear level. On the other hand it was shown 
recently by de Rham and Gobadadze in ^ that it is possible to find such a formulation of 
the massive gravity which is ghost free in the decoupling limit. Then it was shown in || 
that this action that was written in the perturbative form can be resumed into fully non- 
linear actions 2 . The general analysis of the constraints of given theory has been performed 
in [a]. It was argued there that it is possible to perform such a redefinition of the shift 
function so that the resulting theory still contains the Hamiltonian constraint. Then it 
was argued that the presence of this constraint allows to eliminate the scalar mode and 
hence the resulting theory is the ghost free massive gravity. However this analysis was 



questioned in |36| where it was argued that it is possible that this constraint is the second 
class constraint so that the phase space of given theory would be odd dimensional. On 
the other hand in the recent paper very nice analysis of the Hamiltonian formulation 
of the most general gauge fixed non-linear massive gravity actions was performed with an 
important conclusion that the Hamiltonian constraints has zero Poisson brackets. Then 
the requirement of the preservation of this constraint during the time evolution of the 



1 For review, see 0]. 

2 For related works, see @, §, [if, 0, @ |u| g |§ g || || [|§ |r[ @, [§, |f, % g || ||. 
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system implies an additional constraint. As a result given theory has the right number 
of constraints for the construction of non-linear massive gravity without additional scalar 
mode 3 . 

All these results suggest that the gauge fixed form of the non-linear massive gravity 
actions could be ghost free theory. On the other hand the manifest diffeomorphism in- 
variance is lacking and one would like to confirm the same result in the gauge invariant 
formulation of the massive gravity action using the Stiickelberg fields. In fact, it was argued 
in that for some special cases such a theory possesses an additional primary constraint 
whose presence implies such a constraint structure of given theory that could eliminate 
one additional scalar mode. The generalization of given work to the case general metric 
was performed fl3|| . We found the Hamiltonian form of given action and determined pri- 
mary and secondary constraints of given theory. Then we showed that the algebra of the 
secondary constraints is closed and takes the form of the algebra of constraints of diffeo- 
morphism invariant theory. Unfortunately we were not able to find an additional constraint 
that could lead to the elimination of the additional scalar mode with exception of the two 
dimensional case where however two dimensional gravity is trivial. 

These results suggest that it is still very instructive to analyze the non-linear mas- 
sive gravity further in order to fully understand to it. For that reason we perform the 
Hamiltonian analysis of the model of non-linear massive gravity action in the form that 



was presented in 38, 12 1. We consider this action written in manifestly diffeomorphism 
invariant way using the collection of the Stiickelberg fields. We perform the Hamiltonian 
analysis and we identify one important fact related to the symplectic structure of given 
theory. Namely we find that the rank of the matrix of the Poisson brackets between auxil- 
iary fields varies over the phase space. Explicitly, for generic point of the phase space when 
the determinant of the matrix V AB = g l Wi4> A dj(f> B is non-zero we find that the theory pos- 
sesses three spatial diffeomorphism constraints and one Hamiltonian constraint that are the 
first class constraints. Then the constraint structure of given theory at generic point of the 



phase space obeys the basic principles of geometrodynamics [39, 42, 43 1. This result also 



implies that for these generic points there is no additional constraint that could eliminate 



one additional scalar mode and this fact agrees with the results derived in 13]. However an 
interesting situation occurs in the subspace of the phase space where det V AB = 0. In this 
case we find that the Poisson bracket between primary constraints Pab ~ 0, where Pab is 
momentum conjugate to the auxiliary field & AB and corresponding secondary constraint 
^ ab degenerates and possesses eigenvectors whose number is determined by the rank of 
the matrix V AB . Then we perform the Hamiltonian analysis of the system in case when 
det V AB = and we find that now there is one additional first class constraint given as 
the linear combinations of the constraints Pab ~ 0. On the other hand we find that this 
secondary constraint is preserved during the time evolution of the system when the lapse 
function N which is the Lagrange multiplier of the Hamiltonian constraint vanishes. As a 
result the Hamiltonian constraint is the second class constraint. However we should stress 
one important difference between diffeomorphism invariant systems and the system where 



3 For related work, see 
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the Hamiltonian is not given as the linear combination of the constraints. To see this let us 
imagine that the equation that determines the time evolution of given constraint reduces 
to the Poisson bracket between constraint and the Hamiltonian. Then in the case when the 
Hamiltonian does not vanish on the constraint surface the requirement of the preservation 
of given constraint would imply an additional constraint. On the other hand in case of the 
diffeomorphism invariant system we find that the Hamiltonian is given as the linear combi- 
nation of the Hamiltonian constraint and spatial diffeomorphism constraints and hence the 
non-zero Poisson bracket between any constraint and Hamiltonian constraint would imply 
the equation which determines the value of the Lagrange multiplier which is the lapse N in 
this case. In summary we find that the number of the physical degrees of freedom defined 
on the phase space subspace det V AB = is the same as the number of physical degrees of 
freedom determined for the generic point of the phase space. 

At this point we should stress one important point. Systems where the Poisson brackets 
between constraints varies along the phase space are not well understood and it is possible 
that they possess very reach dynamics 4 . It is not completely clear for us how to proceed 
with such systems and for that reason we should be very careful with giving some definitive 
conclusions. It would be certainly very interesting to study the phase space structure of 
given theory in more details, as for example in the paper [46]. Explicitly, it is not clear how 



to proceed with the theory on the degenerative surface det V AB = 0. Shall we interpret 
det V AB = as an additional constraint that eliminate one additional degree of freedom? 
Such a condition can be imposed by the projectable condition <fP = (fP{t) which clearly does 
not correspond to any local degree of freedom. Further, it would be also very interesting 
to analyze the flow of the system near the surface det V AB = 0. We leave all these open 
questions for future. 

The structure of given paper is as follows. In the next section @) we give an overview 
of the non-linear massive gravity in the formulation presented in |38|, 0|. We perform 
the Hamiltonian analysis of given theory and we identify the primary and the secondary 
constraints for the generic point of the phase space where det V AB ^ 0. Then in section 
@ we analyze the theory on the phase space subspace det V = and discuss issues that 
are related to such systems. 

2. Non-Linear Massive Gravity 

Let us consider following non-linear massive gravity action || 

S = Ml f d 4 xy^^R(g) - - A M 2 p m 2 f (fx^gU^f) ■ (2.1) 

Note that by definition g^ v and fu U transform under general diffeomorphism transforma- 
tions x'^ = x ,fJ, (x) as 

dx'V dx' u dr p dr a 



4 For some works, see Q Eg), Egi. 
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Now the requirement that the combination g^ 1 f has to be diffeomorphism invariant im- 
plies that the potential IA has to contain the trace over space-time indices. Further, it 
is convenient to parameterize the tensor fu V using four scalar fields (j) A and some fixed 
auxiliary metric so that 

Uu = d^ A d u 4> B ] AB {4>) , (2.3) 

where the metric Jab is invariant under diffeomorphism transformation x'^ = x fJ '(x') which 
however transforms as a tensor under reparemetrizations of 4> A . In what follows we consider 
Fab = VAB, where tjab = diag(-l, 1, 1, 1). 

The fundamental ingredient of the non-linear massive gravity is the potential term. 
The most general forms of this potential were derived in |[ . Let us consider the minimal 
form of the potential introduced in M 

U(g,H) = -4 (</C> 2 - 




n>l } n>2 



(2.4) 



where we now have 

oo 

K$ = 5% - Vtf -hh = -^2 d n (H n y v , d n 

(2.5) 

and where {H n f v = H^H^ . . . H^ n ~\ 



(2n)! 
(1 -2n)(n!) 2 4 r 



As was observed in [13] in order to find the Hamiltonian formulation of given action it 
is more convenient to consider the potential term written as the trace over Lorentz indices 
rather then the curved space ones. Explicitly, we have 

U = -4(< t > 2 - < t 2 >) = -4(/C A A ) 2 + ±t A B t\ , 

(2.6) 

where we defined 



K% = S A B - ^5 A B - A A B . (2.7) 

Even if it is possible to perform the Hamiltonian analysis for general form of the potential 
we now consider following potential term 

L ma tt = -2M 2 m 2 ^NTry/X . (2.8) 



Following [11, 12] we introduce auxiliary fields B and rewrite the potential term into the 
form 

Cmatt = -M 2 m 2 ^N(^ A B + ($-y B A B A ) , A B A = g^d^ B d v ct> A . (2.9) 
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Note that we have to presume an existence of the inverse matrix (^ 1 ) A B . Then in order 
to see the equivalence between ( |2,8| ) and ( |2.9|) let us perform the variation of (|2.9| ) with 
respect to <& A B 

B 

that implies 



^-YcA^-Yb = (2.10) 



= A A C =► * B = V AA b ■ (2- 11 ) 
Then inserting (2.11) into ( |2.9| ) we obtain ( |2.8| ) and hence an equivalence between these 
two formulations is established. 

Now we are ready to proceed to the Hamiltonian formalism of given theory. Explicitly, 
we use 3 + 1 notation |39| 5 and write the four dimensional metric components as 

goo = -N 2 + Nig %3 Nj , g Qi = N { , g^ = g tj , 



N 2 ' N 2 ' N 2 



(2.12) 



Note also that 4— dimensional scalar curvature has following decomposition 

(4) i? = KijQ m K u + , (2.13) 
where is three-dimensional spatial curvature, K{j is extrinsic curvature defined as 

Kij = ^( d t9ij ~ ViNj - VjNi) , (2.14) 

where Vj is covariant derivative built from the metric components g\j. Note also that Q l i kl 
is de Witt metric defined as 

gijkl = l {gikgj l + g il gjk) _ gl3g kl g.. ki = l( g . kgji + mg , k) _ l gijm (2 15) 

Finally note that in ( 2.13[) we omitted terms proportional to the covariant derivatives which 
induce the boundary terms that vanish for suitable chosen boundary conditions. Then in 
3 + 1 formalism A A B takes the form 

A A B = -V n <p A V n <p B + g^d^d^B , V n A = ^{dt<f> A - N'di^) . (2.16) 



Using this notation we rewrite the mass term ( |2.9[ ) into the form 

Cmatt = -M*m 2 T/gN($ AB VBA - (S^UsV^Vn^ + {$~ 1 )abV ba ) , 
V BA = g^d^d^ , $> AB = VAC 1> C A . 

(2.17) 

The Hamiltonian analysis of the gravity part of the action is well known. Explicitly, the 
momenta conjugate to N, N l are the primary constraints of the theory 

7rjv(x) « , 7Ti(x) « (2.18) 



3 For review, see E(| . 
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while the Hamiltonian takes the form 



n GR = Nn GR + N i n GR ; 
nf R = -2g tk V^ k , 



(2.19) 



where IT 1 ! are momenta conjugate to gij with following non-zero Poisson brackets 

{fly(x),7r«(y)} = \ (8*5) + %S}) <5(x - y) . (2.20) 
Finally note that ttn^i have following Poisson brackets with N, N % 

{AT(x), ^(y)} = <5(x - y) , (iV^x), ^(y)} = 8)5{x - y) . (2.21) 

Now we proceed to the Hamiltonian analysis of the scalar field part of the action. The 
momenta conjugate <& AB are the primary constraints 

Pab = fi^ * ' (2 - 22) 



while the momenta conjugate to <p A are related to the time derivative of <f> A by 

PA = = 2M> 2 ^(<I>- 1 )^V n ^ . (2.23) 

Note that these canonical variables have following non-zero Poisson brackets 



{^ B (x),P CD (y)} = ^cSd + SdSBMx-y) , {^(x),p s (y)} = ^«5(x-y) . (2.24) 
Using these variables we find that the Hamiltonian for the scalar fields takes the form 

n sc = Nuf + ^nf + n AB p BA , uf = PA di(p A , 

= ^^ ^PBPA + Mlrr?^ A * mA + {^) AB V BA \ . 

(2.25) 

Now we analyze the requirement of the preservation of the primary constraints during the 
time evolution of the system. Note that the Hamiltonian is equal to 



H 



= J d 3 x{NH T + WHi + n AB P BA + u N ir N + uVi) , (2.26) 



where Ht = Hj* R + = + ^f c an d where O" 45 ,^ and u l are Lagrange 

multipliers corresponding to the primary constraints P AB ~ 0, 71^ and 7Tj 0. 

As usual the requirement of the preservation of the primary constraints ttn ~ 0, 7Tj « 
implies following secondary constraints 

Ht ~ , » . (2.27) 
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On the other hand the requirement of the preservation of the primary constraints Pab ~ 
implies following secondary ones 

N 

OtPAB = {PAB,H} = - 4M2pm2V gPAPB + 

+ NM^m 2 ^g[-r, B A + {^ 1 )acV CD (^- 1 )db\ = N* AB « 

(2.28) 

using 

{P j4B (x),($- 1 )cD(y)} = -(^ 1 )c^(y){^B(x),^ L (y)} (S-^LuCy) = 

= + 1 )AD(f- 1 )Bc)(x)<5(x-y) . 

(2.29) 

For further purposes it is useful to know the Poisson bracket between Pab and fco 
{P AB (x),^ CD (y)} = ^m 2 ^^ -1 )^* -1 ^ + {^M^bcW^ '(^fd + 

+ (^hEV^i^AF^BD + (^AD^-^BFMx ~ y) = 

= Aab,cd(x)5(x - y) . 

(2.30) 

Note that there are non-trivial Poisson brackets {^ab(x), ^cD{y)} and {^ab(x), ?£r(y)} = 
A j 4b( x , y). Fortunately their explicit forms will not be needed. Finally note that it is useful 
to write the total Hamiltonian in the form 

H T = J d^{NU T + n AB P BA + T AB ^ BA ) + TsiN*) , 

T s (N l ) = J d^N l (n t + d^ AB P BA ) . 

(2.31) 

Using the smeared form of the diffeomorphism constraint Ts(N l ) we find the Poisson 
bracket 

{TsiN*), Z AB } = -diN^AB - N%Z A b ■ (2.32) 

Then it is easy to determine the time evolution of all constraints. Note also that the Poisson 
brackets between smeared form of the Hamiltonian and diffeomorphism constraints take 
the standard form 

{T T (N),T T (M)} = Ts{(NdjM - MdjN)^) , 
{T s (iV i ),T 5 (M- ? ')} = T s (N^d j M i - M'djN*) , 
{T S (N%T T (N)} = T T (d k NN k ) . 

(2.33) 
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In case of the constraints Pab and X P ab we obtain 
d t P AB = {Pab,H t } « A AB ,cdT cd = , 

flt*AB « / d 3 x({* AB ,^ r (x)} + r CD {^ i? ,^ D (x)} + CD {^ i? ,P cl) (x)}) = . 

(2.34) 

The crucial point is to find the solution of the equation 

A AB ,cdT cd = . (2.35) 

To find such a solution we consider an ansatz 

T CD = u c u D , (2.36) 

where u c are 4— dimensional vectors. Then the condition ( |2.35|) implies 

C*- 1 )^* = . (2.37) 

By presumption $ is non-singular matrix so that det(<J> _1 ) ^ 0. Then clearly the only 
solution of the equations ( 2.37| ) is u A = so that T CD = 0. However there is also another 



possibility. In fact, let us consider the ansatz 

T CD = ^CA^DB^ (238) 

so that the condition Aab,cd^ CD = is equivalent to the equation 

V AB u B = . (2.39) 

For generic point of the phase space we can demand that det V AB ^ and hence the above 
condition is solved by ua = 0. However we should stress that the case detV AB = is 
also very important and deservers separate treatment that will be performed in the next 
section. 

Returning to the non-degenerative case we find that thanks to the solution Y AB = 
the second equation in Q2.34 ) takes the form 



J d 3 xiV(x) A AB (x, y) - Q CD (y) A CDiAB (y) = . (2.40) 

Again, due to the properties of the matrix Aab,cd we find that it has unique solution. As 
a result the Lagrange multipliers T CD and £l CD are determined while A r (x) is free which is 
the reflection of the fact that the temporal diffeomorphism is preserved. To see this from 
different point of view note that the constraint T~Lt is preserved during the time evolution 
of the system since 

Wt = {Ut, H t } « J d 3 xr AB (x) {H T , *ab(x)} = (2.41) 



-8- 



due to the fact that T AB = 0. As a result we have ten second class constraints Pab 
0, ^ ab ~ 0. Solving the constraint ^ ab = we find 



= sj ( 4(i ff> + V»(V-') M . (242) 
Then inserting this result into the scalar field part of the Hamiltonian constraint we obtain 



W = 2Mj™V^ ( 4(J ff> + Vac) VC* . (243) 
This form of the Hamiltonian constraint agrees with the Hamiltonian constraint found 



recently in [13]. 



Let us outline our results. We analyzed the model of non-linear massive gravity with 
Stiickelberg fields and with also auxiliary fields that allow to consider theory without the 
square root structure. We shown that the structure of the constrains has the standard 
form corresponding to the diffeomorphism invariant theory without any new additional 
constraint that could eliminate one scalar degree of freedom. 

However the analysis presented in this paper was based on one important presumption 
which is the requirement that det V AB ^ 0. In fact, when we try to impose the unitary 
gauge 

G A : <f> A - x»5 A = (2.44) 

we find that V 00 = V 0i = for i = 1,2,3 and consequently detV AB = 0. In other words 
the unitary gauge ( p. 44 ) cannot be imposed at the generic points of the phase space where 



det V AB 7^ 0. Clearly the analysis of the system where det V AB = deserves separate 
treatment. 



3. Degenerative Case det V = 

As we argued above the non-linear massive theory has an interesting property that its 
phase space structure changes when we move over the phase space. To see this note that 
for the subspace of the phase space characterized by the condition det V AB = we find 
that the matrix Aab,cd possesses the zero modes. Explicitly, for the ansatz 

T AB = ^AC Uc ^BD Ud = T BA ^.l) 

we obtain 

A A b,cdF CD = Mlm 2 ^[u A (^- l )BEV EF u F + (<S>- l ) AE V EF u F u B ] 

(3.2) 

and hence this equation is equal to zero for ua obeying the equation 

V AB u B = . (3.3) 
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Clearly the number of zero modes is determined by the rank of the matrix V AB . We 
restrict ourselves to the case when the number of zero modes is equal to one and we denote 
corresponding vector as w\. Let us then consider following linear combination 

pM = q>AC u c QBD u O D p AB (34) 

Then {P M , Sab} ~ and using ( |3.3| ) we find 

8 t P M = {P M ,H T } »0 (3.5) 

and hence P M is the first class constraint of the theory. In other words from 10 constraints 
P AB we find one the first class constraint and 9 the second class constraints. Note that these 
second class constraints are preserved during the time evolution of the system on condition 
when corresponding Lagrange multipliers are zero. Let us consider the constraint 

Vl/M = ^ab^ AC u° c ^ bd u° c . (3.6) 

Then the time evolution of given constraint is equal to 

d t ^M « {*M,H T } = J ( f 3 x(iV(x){^M,?Mx)} + 

+ {* M ,^ B (x)}r AB (x) + {y M ,p AB (x)}n AB (x)) = 
= J d^(N{^ M ,n T }) = o 

(3.7) 

using the fact that {*m(y), ^ab(*)} T ab {x) = {ty M (y), *m(x)} = and also using the 
fact that 

{*m(x), Pab(y)} ~ u° c $ CE u° D $ DF {^f(x), Pab (y)} = . (3.8) 
To proceed further we have to calculate the Poisson bracket between *m = j^p 
M^m 2 ^grj AB ^ AC ^ BD u° c u D and U T 

{y M (x),n T (y)} = 

' Tgij ^ l (u° A ^ A B PB ) 2 5(x - y) + l - giJ ^\ AB ^ Ac u° c ^ BD u%5^ - y) - 



- (u^pB^g^d^djd^ - y) - -i=( U ° ^B^rti [V&°A-^*Bc] <K X " y) 

(3.9) 

using 

{VsW.^)} = ^^^(x)S(x " y) • (3 - 10) 

Since (3J3) is non-zero we obtain that the equation (3/7) implies that the Lagrange multi- 
plier N should be equal to zero. 

However at this place we should stress one important point. Let us imagine that 
instead of the Hamiltonian constraint we have Hamiltonian that does not vanish on the 
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constraint surface. In this case the equation (|3.7| ) would imply an additional constraint that 
would be the second class constraint together with \& m- However the crucial point is that 
the original Hamiltonian is the linear combination of the constraints and consequently the 
non-zero Poisson bracket ( |3.9| ) determines corresponding Lagrange multiplier rather then 
an additional constraint. As a result ^> m and %t are the second class constraints. We 
also have 18 the second class constraints from original VP ab and Pab- These constraints 
together with the first class constraint Pm eliminates Pab and & AB while Ht and 
eliminates one scalar degree of freedom. The remaining 3 scalar degrees of freedom can be 
eliminated by gauge fixing Hi. As a result we obtain that the physical degrees of freedom 
are 6 metric components gij. 

At this place we should stress that we have to be very careful with such conclusion. 
The problem is that it is not completely clear to us how to proceed with the system that 
degenerates at the subspace of the phase space. More precisely, it is not clear whether 
we should interpret the condition det V AB = as the primary constraint of the theory. 
However the presence of such a constraint would completely change the constraint structure 
of given theory. We leave the complete analysis of this situation for future. 

In any case it is very interesting problem that the Poisson brackets between constraints 
depends on the phase space variables. It is possible that this fact could provide an expla- 
nation why the gauge fixed form of the theory seems to be able to eliminate one additional 



scalar mode [37, while in the case of the manifestly diffeomorphism invariant theory 
defined at the generic point of the phase space such an additional constraint is lacking. 
Acknowledgements: This work was supported by the Grant agency of the Czech republic 
under the grant P201/12/G028. 
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